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1 INTRODUCTION 



ABSTRACT 

We have detected the rare phenomenon of stable, drifting sub-pulse behaviour in two 
pulsars discovered in the recent Swinburne intermediate latitude pulsar survey. The 
pulsars; PSR J1231-47 and PSR J1919+0134, have approximate periods (P) of 1.873 
and 1.6039 seconds respectively. 

Both pulsars have multi-component profiles, and distinct drifting is observed 
across them. We have identified a single drift mode in both pulsars: the drift rate 
for PSR J1231-47 being 5.4(1) ms P'^ and 5.8(2) ms P'^ for PSR 1919-^0134. The 
drifting is linear across the profile with no departure from linearity at the edges within 
the sensitivity of our observations. 



drift rate of an average drift band as a function of pulse 
phase. 



The Swinburne intermediate latitude survey ( Edwards et al 
2001 ) has discovered 69 pulsars, 8 of which are recycled. 
In the process of confirming candidates a number of pul- 
sars were found to display interesting emission behaviour, 
namely pulse nulling and drifting. Two pulsars in particu- 
lar, PSR J1231-47 and PSR J1919-H0134, exhibit regular 
drifting sub-pulses. 

Drifting sub-pulse behaviour has been considered a lit- 
mus test for models of pulsar emission since the discovery 
of such periodicities in pulsar emission by Drake and Craft 
(1968). Although many pulsars display sub-pulse intensity 
variations few present orderly persistent drifting sub-pulses. 
It is considered by Rankin (1986) to be a purely geometrical 
effect. Assuming the intensity variations are due to conal 
sub-beams, then drifting will be apparent if the pulse pro- 
file is a result of an almost tangential cut of the line of 
sight across the cone. This implies that drifting sub-pulses 
would be more apparent in conal single profiles and that well 
resolved double profiles should display intensity variations 
that do not drift. The newly discovered pulsars are both 
members of an interm ediate profile class, that of "barely re- 
solved conal double" (Rankin 1983). They also display the 
very rare property of a stable drift pattern. 

We present both 2 dimensional auto-correlation anal- 
ysis (Vivekanand and Joshi 1997) and fluctuation spectra 
(Backer 1973; Taylor Manchester and Huguenin 1975) for 
both pulsars; providing an initial characterisation of their 
sub-pulse behaviour. We also present an examination of the 



2 OBSERVATIONS 

All the observations presented here were taken with the 
Parkes 64 metre radio telescope using the central beam of 
the 13 beam multi-beam receiver and the 96 x 3 MHz filter 
bank (Staveley-Smith et al 1996), at a central frequency of 
1374 MHz. Both polarisations were summed together and 
the result one bit sampled. The sample rate varies between 
the observations, and is 125 fis for PSR J1919-I-0134 and 
500 for PSR J1231-47. This time series is folded into 512 
phase bins at the topocentric pulsar spin period. 

Each period can be represented by a row in a 



"longitude-time" plot ( 


Taylor, Manchester & Huguenin 


1975; 


Manchester & Taylor 1977 


). Two such plots are shown 



These arrays form the basis of the following analysis. 

Many different methods have been employed to exam- 



ine drifting behaviour, including cross-correlation ( Proszyn 
ski & Wolszczan 



1986) 



Fourier phase methods and fluctu- 
ation spectra (Backer 1973, D eshpande and Rankin 2001 ) 
and auto-correlation methods (Vivekanand & Joshi 1997). 



Both fluctuation spectra and auto-correlation analyses are 
presented here. 
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Figure 1. Longitude— time plot showing the drifting sub-pulses of J1231+47 and J1919+0134. The contour levels have been set to a 
scale that increases the contrast. Phase has been set arbitrarily. The dotted line indicates the level of spectral power in the fluctuation 
spectrum at that longitude. 



Table 1. Astrometric, spin and dispersion parameters for the pulsars under investigation, from Edwards ct. al. (2001). It should be 
noted that there is no complete timing solution for PSR J1231-47. 



Name 


R.A. (J2000) 


Dec. (J2000) 


P 


P Epoch 


P 


DM 








(s) 


(MJD) 


(10-1^) 


(cm"^ pc) 


J1231-47 


12'^31™40(50) 


-47° 46(7)' 


1.87304(2) 


52056.3 




31(3) 


J1919+0134 


19'^19™43!62(3) 


+01°34'56'.'5(7) 


1.60398355528(6) 


51650.0 


0.589(11) 


191.9(4) 



2.1 Fluctuation spectra 

Fluctuation spectra were formed by performing a spectral 
analysis of t he lon gitude-time data. The process is discussed 
in Backer (197J). The Fourier transform was performed 
upon each longitudinal column to produce an indication of 
the level of periodicity as a function of pulse phase. The 
spectra were then normalised by first subtracting the power 
spectra of an off-pulse region from all the on-pulse spectra. 
Each spectrum was then divided by the square of the aver- 
age signal level in the longitudinal column from which the 
spectrum was formed. Figure |^ shows the normalised aver- 
age pulse profile along with the normalised spectral power 
as a function of longitude measured at the peak of the fiuc- 
tuation spectrum. Both pulsars are in fact very similar. The 
fiuctuation, although great in longitudinal extent, is concen- 
trated in the stronger of the two components. 

Analysis of the spectra themselves provided values for 
the frequency of the amplitude variations. The values ob- 
tained for the frequency of the main spectral feature were: 
0.050±0.002 (cycles/period) for PSR J1231-47, and 0.156± 
0.008 (cycles/period), for PSR J1919+0134. 



In order to further investigate this phenomenon, a 1 
hour integration of PSR J1231-47 has been obtained. Due 
to interstellar scintillation effects the 1 hour observation had 
a signal to noise ratio such that sub-pulses are not directly 
observable in a longitude-time plot. Furthermore it suffered 
from a burst of interference approximately half way through. 
In order to mitigate the interference effects the observation 
was broken into two sub-integrations. Nevertheless fluctua- 
tion spectra formed from these sub-integrations displayed a 
significant spectral feature at 0.0528±0.0008 (cycles/period) 
and 0.0515 ± 0.0004 (cycles/period) respectively. This fea- 
ture is consistent with that found in the shorter observation. 
The quoted errors are under-estimates being simply half a 
frequency bin- width. These values are consistent with those 
presented in Table bl 



2.2 Average properties 

In order to examine the features of a "typical" drift band 
some sort of band averaging had to be performed. This was 
achieved by the construction of a two dimensional auto- 
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correlation function (2DAC F) of the longitude - time plots 
(Vivekanand & Joshi 1997). The auto-correlation function 
was formed in a number of stages. A region of the longi- 
tude time frame was chosen. The region is composed of the 
on-pulse region in longitude domain and the length of the 
observation in the time domain. The mean value of this sub- 
array was then subtracted from each point. The purpose of 
this was to reduce the height of the zero-lag spike in the 
auto-correlation. The array was then padded on all 4 sides 
with a number of zeroes commensurate with the number of 
lags required in the 2DACF. A forward Fourier transform 
was then performed. The resultant was then multiplied by 
its complex conjugate. The product was then inverse trans- 
formed, producing the 2DACF. 

The 2DACF provides a method for examining the "av- 
erage" properties of the drifting bands by providing the cor- 
relation of each drift band with itself and its neighbours. We 
can examine the evolution of the drift rate by comparing the 
slopes and separations of the correlated bands. The values 
for P2 (the p hase separation), and P3 (period separation) 



(Backer 1973) can be found by analysis of straight line fits 
to the peaks in the lag-lag plane. The straight line fits were 
performed on the central regions of each drift band, where 
the signal to noise ratio was highest. In order to obtain the 
parameters of the straight line fit we assumed the data was 
adequately modeled by a straight line and an initial fit was 
performed. Then measurement errors of magnitude 1 stan- 
dard deviation were ascribed to each point, the fit performed 
again and the chi-square statistic minimised to obtain val- 
ues and errors for the fitted parameters. This method is 
described in Press et al. (1986). 

The values for the separation between pairs of bands in 
time (P3) and phase (P2) are given in Table ^. The values 
were obtained by combining the values for individual drift 
bands in the lag-lag plane obtained from the straight line 
fits. 



2.3 Drift rate across the profiles 

Many pulsars which display the property of drifting sub- 
pulses exhibit a change in sub-pulse drift rate across the 
pulse pr ofile. This property i nay be used to test e mission 
models ([Krishnamohan 1980|; |Oster fc Sieber 1976|; [Wright 



1981) 



We have examined the straight line fits to the peaks in 
the lag-lag plane in an attempt to identify any clear depar- 
ture from linearity. No discernible, regular structure indica- 
tive of a changing drift rate is apparent in the residuals. Al- 
though it should be noted that the "grand averaging" prop- 
erties of the 2DACF method will remove any non-linearity if 
the drift rate varies between bands. The low signal to noise 
ratio of the sub-pulses at the edges of the drift band make 
investigation of the band behaviour in this region difficult. 



Table 2. The parameters of a "representative" or "average" 
drift band for the two pulsars. The parameters were derived from 
analysis of the straight line fits to peaks in the 2 dimensional 
auto-correlation functions. P2 is the separation of two drift bands 
in phase. P3 is the separation in periods of the drift bands. The 
drift rate is given as the number of periods required for an average 
drift band to cross unit phase. The drift rate in milliseconds per 
period is also given. 



Name 

J1231-47 

J1919-I-0134 



P2 (</-) 

0.0557(9) 

0.023(1) 



P3 (P) 
19.3(2) 
6.5(2) 



Drift rate (P/(/>) 

348(8) 

277(8) 



(ms/P) 

5.4(1) 

5.8(2) 



gitudinal extent of the drift is interesting; it is rare for pul- 
sars of this profile class to display a measurable drift acros s 
the whole profile ( Rankin 1986 ; Hankins fc Wolszczan 1987 ). 
Both of these pulsars display a measurable P3 across the ma- 
jority of the profile and the drift is continuous. Furthermore 
it is rare for pulsars of this profile class to show such a clear 
sense of drift, as simple periodic intensity variations at fixed 
longitudinal positions are more common. 

It appears that these pulsars are rare and exemplary 
exponents of the drifting sub pulse phenomena; as such these 
pulsars will be of great benefit in investigations of the radio 
emission mechanism in pulsars. 
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3 DISCUSSION 



Very regular, stable, persistent drifting pulsars are rare. As 
such these pulsars represent a welcome addition to the pulsar 
family. The initial observations are intriguing, the linearity 
of the drift ma y provide some information regar ding emis- 
sion geometry ( Krishnamohan 1980 ; Wright 1981 ) . The lon- 
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Figure 2. The 2 dimensional autocorrelations of the longitude-time information. A straight line has been fit to the central peaks in the 
lag-lag plane in order to characterise the drift band behaviour. The parameters are given in Table Gl 
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